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Abstract — In this paper , We have extended Cohen’s 

work[l] to approximate cos( 1 /2 nx4) and sin(‘/2 nx4) , 

0 < x < 1 , in terms of Chebyshev polynomials . We have used 
two approaches to effect the agreement in the result of 
computations for Chebyshev coefficients to 8 digits of accuracy. 
The first approach is the fundamental Chebyshev IDENTITY 
derived by Cohen [1], However, We have DERIVED this 
IDENTITY again in THEOREM 1 and We have observed that 
each entry element in Table 2 DERIVED by Cohen’s work [1] 
and listed in Appendix A can be evaluated separately. 

This , in contrast to COHEN’S DERIVATION [1] in 
which each individual entry in Table 2 can not be evaluated 
separately unless the evaluations of the previous rows are 
calculated first . The second method, we have apply is the direct 
series summation techniques. 

Index Terms — Chebyshev Identities, Clenshaw 


I. INTRODUCTION 

The BACKGROUND of this paper is that Clenshaw [4] has 
approximated 

sin Q/i nx) & cos Q/i nx) , -1 < x < 1 as an expansion of 
Chebyshev polynomials . The coefficients are determined 
exactly in terms of Bessel functions and then approximated 
numerically to 20 digits of accuracy as in [4], Cohen [1] has 
extended Clenshaw’ s work [4] to approximate cos (V 2 nx 2 ) & 
sin (V 2 nx 2 ), 

0 < x < 1 as an expansion of Chebyshev series. The 
coefficients are determined to 20 digits of accuracy using two 
approaches. 

The first approach he has derived a 
fundamental Chebyshev IDENTITY for expressing 
Chebyshev polynomials T r (x 2 ) as a linear approximation in 
Chebyshev polynomials Tj(x), 

r, j=0(l) n. This fundamental IDENTITY that We shall 
derived again produces Table 2 in Appendix A. The second 
approach in [1] was the continued fraction technique. One 
APPLICATION of this, demonstrated, 
in Cohen's work[l] is the evaluation of FRESNEL 
INTEGRAL . 

In Cohen’s work [1], Chebyshev polynomials approximations 
are expressed as 

cos V2 nx 2 = XB 2r T 2r (x) (1) 

and 

sin 1 /2nx 2 = X(-l) r 8 2r T 2r (x) (2) 

where the coefficients B's are given in Table 1. The 
evaluations of the coefficients B's in Table 1 are determined 
by establishing Chebyshev identities for T r (x 2 ) as described 
by Cohen [1], 
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That is 

r 

2 r T r (x 2 ) = £ a r , k T 2r _ 2k (x) (3) 

k = 0 

where the coefficients a’s are listed in Table 2. 

We observed that the determination of a single entry in 
Table 2 is not feasible in Cohen’s derivation [1] unless the 
evaluations for the entries start from the first few rows. We 
overcome such a difficulty by 

Table 1 Coefficients in the expansion for sinVinx 2 
and cosVinx 2 (0< x < 1) 


cos V2TIX 2 = X B 2r T 2r (x) , sinViHx 2 = I (-l) r l5 2r T 2r (x) 


2r 15 2r 


0 0.60219 47012 55546 40329 

2 -0.51362 51666 79107 02511 

4 -0.10354 63442 62963 75381 

6 +0.01373 20342 34358 55321 

8 +0.00135 86698 38090 36178 

10 -0.00010 72630 94406 00221 Checks are 

12 - 70462 96793 46857 provided by 

14 + 3963 90250 61486 X(-l) r f5 2r =l 

16 + 194 99597 75588 &Ifl 2r =0 

18 - 8 52292 89262 

20 - 33516 50652 

22 + 1197 93739 

24 + 39 24123 

26 - 1 18639 

28 - 3330 

30 + 87 

32 + 2 


deriving a formula for Table 2 again and thus each entry 
element in it can be evaluated directly. 

H. CHEBYSHEV IDENTITIES 
The first well known identity that we have rewritten , is 
the explicit representation of x n in terms of Chebyshev 
polynomial Tj(x) , j=0(l)n , for some positive integer n. 


That is 


x" = £ cj Tj(x) 

(4) 

J=0 


where Cj = 2' n+1 C(j+n)/2 

(5) 


for j =l(l)n if n is odd & j =0(2)n if n is even . The prime 
dash in the summation means that coefficient of T 0 (x) in (4) 
should be halved . 

The next identity that we have derived again without proof 
is the expansion of Chebyshev polynomials T n (x) as 
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successive powers in x , where n is a positive integer T n (x) 

r 

= £ a j x "~ 2j (6) 

7=0 

where the coefficients a's are determined by 


c%r 


do - 2 ” 1 

and o j = (-l) J 2 n_2j l C”r/ _1 - (7) 

J 7 

for j = 1 (1) r and r = [n/2] is the integer part of n/2 Where 
' n j ^ is the Corresponding Binomial coefficient. For 

an alternative expression of Eqn. (6), see [2], 

The third identity that we shall derive again, is the 
expansion of T n (x 2 ) in terms of Chebyshev series, as shown in 
Theorem 1. For an alternative derivation, see [1]. 

Theorem 1. 

If T n (x) is the Chebyshev polynomial of degree n, 
where n is a positive integer , then we can expand T n (x 2 ) in 
terms of a Chebyshev series as 


T n (x 2 ) = £ a n , k T 2k (x) (8) 

;= o 

where 

»■*- e» 

for k=0(l)n , r = [n/2] and a. t for j = 0(l)r are described by 
Eqn .(7). The prime dash on the sign sigma means the 
coefficients a n>0 are to be halved , for each integer n>l 

Proof 

Eqn. (6) implies that 


T n (x 2 )= X OjX 

7=0 

while Eqn. (5) implies 


( 10 ) 


n-2j 

x 2n-4j = £ Ck T 2k (x 

k = 0 

hence , We can write 


Tn(x 2 ) = £ a n , k T 2k (x) 

k = 0 

with 

2 jr(x 2 )r 2 *«, 

" ,k ' n_ - 1 


r 2 \ x 2H 4j TikW 
a„, k = £ Oj(- J r 2 d 
7=0 n - 1 V 1 ” X 


The method describes the evaluations of the integral 
enclosed within the parentheses can be found in Froberg's 
book[3] , page(79). 

Hence , We get 


a n ,k 


= a. 




2n—4j 

n-2j+k 


(14) 


for k=0(l)n , where 

aj for j=0(l)n are given by Eqn .(7) and r = [n/2]. 

We remark that the elements a n k , for k =0(l)n , are the 
entries for Table 2 derived earlier by Cohen’s 
work[l] . If initial coefficients a n ,o are halved , then , we can 
compute each entry element that lies anywhere inside Table 2 
. For example , further elements in Table 2 are given by 
a 25,23 = 1 125 and a 29 , 27 =1537 


III. APPROXIMATION OF COS ViUX 4 AND SIN ViTlX 4 
We shall first apply chebyshev expansion series for cos 
Vinx 4 , where the general coefficient is evaluated using direct 
series summation technique . That is when 

cos'Anx 4 = £ B r T 2r (x) (15) 

r = 0 

then B,= — [ C0S — — dx (16) 

* i aF J 1 

If We expand cosVinx 4 as 
oo 

cosVinx 4 = ^ (-l) i ( 1 / 2 n) 2j x 8 V( 2 j!) (17) 

7=0 

and substitute x = cos0 , then Eqn. (16) can be rewritten as 

B s obW t!_, 

' M> ,2 J Vl 

2 g j 

* ( — j (cos 6) J cos 2r0 d0) (18) 

7T Q 

hence 

B r = “ (-iy(/2n) 2j 2- 8 J +1 cl J j+ r /m (19) 

7=0 
for r = 0(l)n 

We observe that the terms in the series (19) are decreasing 
rapidly to zero and all the terms are identically zero whenever 
r >4j. 

Alternatively, We knew that if we replace x by x 2 in 
Eqn. (1), We get 

cos Vinx 4 = X B 2r T 2r (x 2 ) (20) 

By using the values of B 2r , r = 0(1)16 , given in Table 1 and 
Chebyshev identities for T r (x 2 ) , We get 

cos 1 /2nx 4 = XB r T 2r (x) (21) 

where the coefficients B's in Eqn. (21) are exactly given by 
Eqn. (19) . The values of the coefficients B's are shown in 
Table 3 rounded to 8 significant digits according to the 
agreement of the evaluation by both methods. 

Similarly, the approximation for sin Vinx 4 can be evaluated 
by direct series summation or Chebyshev identities for T r (x 2 ) 

sin 1 /2nx 4 =XB r T 2r (x) (22) 

where 
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6 r = x (-lye/tfo 2 * 1 ^- 3 c^Jjtt+2 /(2j+ 1} ! (23) 

for r = 0 (l)n 

The values for the coefficients B's given by Eqn.(23) are 
listed in Table 4 rounded to 8 significant digits according to 
the agreement in the results by both methods 

Table 3 The coefficients B’s in the expansion 
for cos Vinx 4 

cos ‘/2nx4 = X I.irT2r(x) 


0 

0.70923853 

1 

-0.45715459 

2 

-0.21256022 

3 

-0.04639916 

4 

+0.00271565 

5 

+0.00352780 Checks are 

6 

+0.00062078 provided by 

7 

+0.00002962 X(-l)rBr=l 

8 

-0.00001443 

9 

-0.00000370 

10 

-0.00000034 

11 

-0.00000001 

Table 4 The coefficients B's in the expansion 


for sin Vinx4 

sin 1 /2nx4 = V BrT2r(x) 


0 

0.45669090 


1 

+0.55985521 


2 

+0.05834041 


3 

-0.05833114 


4 

-0.01536607 

Checks are 

5 

-0.00165872 

provided by 

6 

+0.00031830 

I(-l)r Br = 0 

7 

+0.00013469 


8 

+0.00001688 


9 

+0.00000003 


10 

-0.00000043 


11 

-0.00000008 



IV. COMPUTATIONAL REMARKS 
The accumulations of rounding errors due to the 
limitedness machine's accuracy of I.B.M. Micro computers 
applied with Basic Language , effect the evaluations of the 
coefficients in Table 3 &Table 4 to be taken to 8 significant 
digits . 

However , both the methods described can retain more 
accuracy for the evaluations of the coefficients in Table 3 & 4 
, whenever recomputation is taken on sufficient machine's 
accuracy .The continued fraction technique can be used 
instead of direct series technique . 
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(r) . 
s OLr—k ' 


i expansion of 2Tr(s?) 


12 13 14 15 16 


8 35 -80 

9 70 -126 

10 -126 280 

11 -252 462 

12 462 -1008 

13 924 -1716 
14-1716 3696 

15 -3432 6435 

16 6435 -13728 


72 


48 

400 


-216 216 


312 


117 


-176 605 500 150 20 


-308 -22 528 -891 286 1287 748 187 22 1 

1224 -1424 1263 -264 -1452 1608 2370 1064 228 24 1 

1248 -273 -1274 2847 -2704 -1014 4420 3978 1456 273 26 1 

-4389 5124 -5026 2772 2233 -6744 2324 9576 6251 1932 322 28 1 

-4950 1925 2940 -8637 11310 -4875 -10200 11815 18186 9345 2500 375 30 1 

16016 -18656 19272 -14112 176 18336 -22756 -6688 32528 31648 13432 3168 432 32 1 
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